Topological invariants for spin-orbit coupled superconductor nanowires 
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We show that a spin-orbit coupled semiconductor nanowire with Zeeman splitting and s-wave 
superconductivity is in symmetry class BDI (and not D as is commonly thought) of the topological 
classification of band Hamiltonians. The class BDI allows for an integer Z topological invariant 
equal to the number of Majorana fermion (MF) modes at each end of the quantum wire protected 
by the chirality symmetry (reality of the Hamiltonian). Thus it is possible for this system (and all 
other d = 1 models related to it by symmetry) to have an arbitrary integer number, not just or 1 
as is commonly assumed, of MFs localized at each end of the wire. The integer counting the number 
of MFs at each end reduces to or 1, and the class BDI reduces to D, in the presence of terms in 
the Hamiltonian that break the chirality symmetry. 
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Introduction: Rashba spin-orbit (SO) coupled semi- 
conductors in dimensions c? = 2, 1 with a Zeeman field 
and proximity-induced s-wave superconductivity have re- 
cently attracted a lot of attention [ll-[l7| . Under suitable 
external conditions these systems can support MF exci- 
tations (defined by second quantized operators 7^ = 7) 
whose statistics is non-Abelian. In = 2 the particle-hole 
(p-h) symmetric Bogoliubov-de Gennes (BdG) Hamilto- 
nian of the Rashba-coupled semiconductor is in the topo- 
logical class D 3) 1^ with an integer Z topological in- 
variant which counts the number of gapless chiral Majo- 
rana modes on the boundary. From dimensional reduc- 
tion (i.e., by putting one of the wave- vectors to zero (20|). 
the gapless boundary Majorana modes in d = 2 reduce 
to zero-energy end Majorana modes in a d = 1 nanowire. 
The dimensional reduction argument therefore suggests 
that the number of possible end Majorana modes in a 
SO coupled nanowire should also in principle be an inte- 
ger. Based on this (and more rigorous arguments below) 
we argue that the Hamiltonian of the system in d = 1 
is in the topological class BDI (in contrast to its being 
in class D in d = 2), characterized by a Z topological 
invariant which counts the number of possible end MFs 
in the nanowire. By rigorous arguments we clarify why 
the d = 1 system is in class BDI and has a Z topological 
invariant, construct an algebraic form of the invariant, 
discuss its relation with the more-frequently-used Z2 in- 
variant, and illustrate these points by constructing ex- 
plicit examples. The work here significantly clarifies the 
topological properties of the SO coupled nanowires and 
related systems which have recently become an important 
focus of attention both theoretically and experimentally. 

The topological class of the SO coupled semiconductor 
is analogous to that of a spinless Px + iPy superconductor. 
In d = 2 the spinless Px + ipy superconductor, with bro- 
ken time reversal (TR) invariance (due to the presence of 
i in the order parameter) , is in class D characterized by a 
Z invariant counting the number of gapless chiral Majo- 
rana modes on the boundary 22|,|2J|. It is also possible to 



define a Z2 invariant which only counts the parity of the 
number of boundary Majorana modes 26]. Dimensional 
reduction arguments suggest that the number of possible 
end MFs in a d = 1 spinless superconductor should also 
be an integer and this has recently been shown explicitly 
21J . Therefore, the Hamiltonian should be in the topo- 



logical class BDI with a Z invariant in d = 1. Note, how- 
ever, that d = 1 Hamiltonians in class BDI are supposed 
to be TR-invariant while the spinless Px + ipy supercon- 
ductor explicitly breaks TR symmetry in d = 2. The key 
to this difference is that, in d = 1, the Hamiltonian can 
be made completely real 2l| while it is necessarily com- 
plex in d = 2. Redefining the time-reversal operator only 
in terms of the complex conjugation operator /C, it fol- 
lows that in d = 2 this symmetry is broken (class D) but 
it remains intact in d = 1 (class BDI). More generally, 
the emergence of the reality condition in d = 1 (reality 
of HndG in the present case implies the chiral symmetry 
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given hy S — IC ■ A where A is the p-h transfor- 



mation operator) changes the symmetry class of both the 
spinless p-wave superconductor and the Rashba-coupled 
system from D in d = 2 to BDI in d = 1. 

Despite the fact that HBdG for a d = 1 Rashba-coupled 
nanowire can be made purely real, the Z invariant is 
not computable as a winding number of the Anderson 
pseudo-spin d-vector defining the BdG Hamiltonian as 
HbiIgO^) = d(k).f |22l424| . Here, r is a vector of Pauli 
matrices defined in the p-h space. This is because, in con- 
trast to a spinless p-wave superconductor in d = 1, the 
components of the d-vector for the nanowire are them- 
selves matrices. More generally, Hbclg for a topologi- 
cal superconducting (TS) system in d = 1 can be real 
(thus preserving the chiral symmetry) but can be a large 
27V X 27V square matrix so the components of the d-vector 
are themselves N x N matrices. We show below how to 
compute the integer topological invariant for this prob- 
lem in terms of a generalized pseudo-spin vector and its 
winding number and connect the value of this integer to 
the number of independent MF modes at each end of 
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the nanowire. Note that multiple Majorana modes at 
each end are protected by the chiral symmetry. We fur- 
ther show that the role of the usual Z2 Pfaffian invariant 
25I 26| reduces to determining the parity of the Z invari- 
ant. In the presence of chiral symmetry breaking terms 
(i.e., terms which introduce complex entries in HbcIg) 
the topological class of the system reduces to D even in 
d = 1 and an odd (even) number of MFs at each end 
becomes equivalent to just one (zero) . 

Z invariant for real BdG Hamiltonians in d = 1; 
To understand the difference between complex and real 
Hamiltonians let us start from the Hamiltonian of a spin- 
less Px + ipy superconductor in 0? = 2, 

iJl(k) = (ek - li)Tz + AxkxTx - AykyTy, (1) 

where k is a two-dimensional wave- vector, fi is the chem- 
ical potential, and , Ay are superconducting pair po- 
tentials along the x, y directions, respectively. Here we 
have used the p-h basis (cj^., c_k) and its hermitian con- 
jugate, and the r matrices in Eq. [1] are defined in this 
basis. Writing this Hamiltonian in terms of the Ander- 
son pseudo-spin vector [23] d(k) as -ffi(k) ~ d{'k).f, we 
see that for spinless Px + ipy superconductor in d = 2 all 
three components of d are non-zero. The group of topo- 
logical invariant is then Z which is the relevant homotopy 
group 712(5*^) of the mapping from the two-dimensional 
k space to the 2-sphere of the 3-component unit vector 
d = d/\d\ [l^,!!^. On the other hand, in d = 1, since the 
corresponding Hamiltonian can be made purely real (A^; 
drops out from Eq. ([ij for the system along the y-axis), 
the vector d has only two components. Noting that the 
fc-space now is also one-dimensional, the topological in- 
variant must again be in Z (class BDI) since 7ri(5^) = Z. 
This invariant is simply the winding number. 
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where 9{k) is the angle the unit vector d makes with, 
say, the z-axis on the y — z plane. The winding num- 
ber counts the number of times the 2-component vec- 
tor d makes a complete cycle in its plane as k varies in 
the one-dimensional Brillouin zone. It is clear that only 
with the breakdown of the reality condition of the BdG 
Hamiltonian (i.e., chiral symmetry given hy S = K, ■ K) 
the symmetry class of the spinless p-wave superconduc- 
tor can change from BDI to D (which is characterized by 
a Z2 invariant) even in d = 1. 

In d = 2, 1 the 4 X 4 BdG Hamihonian iJ2(k) of 
a Rashba-coupled semiconductor with Zeeman coupling 
and proximity induced s-wave superconductivity is given 
by, 

ff2(k) = (ek - y)T^ + VzS ■ ctt^ + akxfJyT^ 

- akyffx + AoUyTy, (S) 



where we have used the 4-component p-h spinor 
{u^{r),Uj^(r),v^(r),vi{r)) (with quasiparticle operators 
given by d'^ = I]o-("cr(^)4(^) + VoCair))), and the PauH 
matrices crx^y^zTTx,y,z act on the spin and particle-hole 
spaces, respectively. In Eq. JS]), the vector S is a suit- 
ably chosen direction of the applied Zeeman spin split- 
ting Vz (e.g., 5' = z in d = 2 for k = {kx,ky), S = x 
in d = 1 for k = fc^), /i is the chemical potential, a is 
the Rashba SO coupling constant, and Ag is an s-wave 
superconducting pair-potential. It is clear that in d = 2 
it is not possible to make the Hamiltonian real because 
of the complex Rashba term. In contrast, in d = 1 H2 
can be made purely real (Rashba term can be made real), 
and one can define a pseudo TR operator in terms of K. 
alone. Then, in d = 1 H2 preserves both p-h as well as 
the new 'time reversal' symmetry and hence is in class 
BDI characterized by a Z invariant. Note, however, that 
in contrast to the case of a spinless p-wave superconduc- 
tor, the components of the d- vector in the present 4x4 
Hamiltonian are themselves 2x2 matrices. It is not clear 
how to define the winding number of the d vector or the 
Z invariant for such BdG Hamiltonians which are larger 
than simple 2x2 matrices as in spinless p-wave super- 
conductor. More generally, the BdG Hamiltonian of a 
TS system in d = 1, despite being real (thus preserving 
the chiral symmetry), can be a large 2N x 2N square 
matrix so that the components of the d- vector are N x N 
matrices. We now show below how to compute the inte- 
ger topological invariant for this problem by generalizing 
the concept of the d- vector winding number for arbitrary 
dimensional matrices. 

For real BdG Hamiltonians in d = 1 it is possible to 
construct an integer topological invariant in terms of a 
generalized winding number. To see this, consider an ar- 
bitrary second quantized BdG Hamiltonian (real or com- 
plex) which can be written as |27| . 



H = 



Ai 



(4) 



To write the BdG matrix in Eq. (|4]) we have used the 
p-h basis where = (4,c,), = (cj,^)^, 

i,j are collective coordinates representing both spatial 
coordinates r, r' and spin indices '1, 4-, and summations 
over i and j are implied. In this form, the gap matrix 
Aij is necessarily an anti-symmetric matrix. If the BdG 
matrix is purely real, Hq^^^ — i?o,ij 



^ - and At = A^. = 



-Ai,-. In this case, we can rewrite the BdG Hamiltonian 



in the form (omitting the i,j indices) 



H 



Ho A 

-A -Ho 



(5) 



where Ho is real-symmetric and A is real anti-symmetric. 
Such a Hamiltonian H has the chiral symmetry [l8l - [20j 
defined as 5 = /C • A (with A = Tx ■ IC in this basis). 
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under which the Hamiltonian is invariant. Since in the p- 
h space the matrix H can be written as H = HoTz + iAry 
it can be made purely off-diagonal by a rotation in the 
p-h space by the unitary transformation U = e^^^'^y . It 
follows that the rotated Hamiltonian 



A 




(6) 



is a symmetric Hamiltonian with the matrix A = Hq + A 
being real. 

Fourier transforming to momentum space A{k) satis- 
fies the constraints A{k) — A*{~k), so that A{K) is real 
at the p-h symmetric points k = K, which are given by 
K — 0, ±7r in d = 1. The Hamiltonian in the k-space can 
be written as, 



UH{k)U'^ 



A{k) \ 
A^(-fc) ) 



(7) 



Notice that the existence of a zero eigenvalue of the ma- 
trix A(k) necessarily implies the existence of a zero eigen- 
value of the BdG Hamiltonian. Therefore, as long as all 
the eigenvalues of A{k) are gapped (which implies that 
the determinant of A(k) is gapped), the Hamiltonian is 
also gapped. It follows that a non-zero Det{A{k)) for all 
values of k in the d — \ Brillouin zone indicates a fully 
gapped Hamiltonian. 

Since A{k) is a purely real matrix at the p-h in- 
variant points k — 0,±7r, Det{A{k)), which is in gen- 
eral a complex number for a general value of fc, is also 
purely real for k — 0,±7r. We now write Det{A{k)) 
as Det{A{k)) = \Det{A{k))\ex:p{ie{k)). Note that e(k) 
must be or ±7r at fc = 0, ±7r. Defining the variable 
z{k) =exp(i6'(fc)) = Det{A{k))/\Det{A{k))\, we caTVTvow 
write the expression for a winding number, 



W 



k=0 



^ dz{k) 
'4k)-' 



(8) 



which can only be an integer (i.e. W £ Z) including zero. 
Similar topological invariants for BDI Hamiltonians in 
d = 1 have also been constructed earlier [2^ [29[ This 
integer is a topological invariant because it is not related 
to any symmetry breaking and yet can only change when 
Det{A{k)) goes through zero somewhere between fc = 
and k = TT which indicates a closure of the spectral gap 
and the associated quantum phase transition. 

We now do a consistency check on the topological in- 
variant W. By straightforward algebra we can show that. 



W = 



-i dz{k) 



dk 



9fe[logz(/c)] 



dk. 



27r 
dk 



-Tr 



( A{k)dkA-\k) 



\ A^{~k)dkA^--\~k) 

/ ^TT[TzH{k)dkH{k)-^]^Wu (9) 



where Wi is the topological invariant written for chi- 
ral systems in terms of the zero-frequency single-particle 
Green's functions G{k) — H^^{k) in Ref. Isol. It can also 
be shown 30] that the difference of the integers Wi be- 
tween two gapped topological systems gives the number 
of zero energy modes at a boundary separating them. It 
follows that W gives the number of MF modes at each 
end of a semiconductor nanowire since the end separates 
the chiral nanowire (with integer invariant W) from vac- 
uum (with W — 0). 

Pfaffian Z2 invariant as parity of the Z invariant: 
Now we derive a formula connecting the Z invariant W 
and the Pfaffian Z2 invariant [2^ [2^ more frequently 
used for a semiconductor nanowire. For this, consider a 
BdG matrix HbcIG with a particle-hole symmetry of the 
form TxHbcIG = ~Hgj^QTx (note that A — JC ■ t^) where 
Tx = rj is the symmetric particle-hole transformation 
matrix satisfying t^t* — 1. Then the matrix HbcIcTx is 
anti-symmetric i.e 



[HBdGTx)'^ — TxHg^Q — —HBdGTx 



(10) 



This allows us to define a Pfaffian Pf{HBdGTx) associ- 
ated with the BdG Hamiltonian. 
Now note that 



A{k) 
-A^i-k) 



Det{A{k)) 



(11) 



at A: = 0, TT, where we have used the fact that Det{U) = 
1. Then, the Pfaffian topological invariant of the BdG 
Hamiltonian [261] , which is the algebraic sign of the quan- 
tity Q, 



( A{k = 0) 

^^\~A^{k = Q) 



A{k ^ tt) 

-A^{k = tt) 



(12) 



is clearly equal to the sign of the Det{A{k = 
0))/Det{A{k = tt)). The sign of Det(A(fc = 
0))/ Det{A{k — tt)) is in turn equal to the parity of W, 
because, from Eq. \8\ 



sign 



Det{A(k = 7r)) z{k = tt) 
Det{A{k = 0))^ ~ z{k = 0) 



(-1) 



w 



(13) 

It follows that the familiar Z2 Pfaffian invariant of the 
d = 1 systems is simply the parity of the more general Z 
invariant of a chiral Hamiltonian. 

Z invariant for spinless superconductors in d = 1: 
Next we show that W as defined in Eq. ([8]) reduces to 
the invariant N (Eq. [2|) for a c? = 1 spinless p-wave super- 
conductor. In this case, the BdG Hamiltonian is given 
by, HBdG = (efc - IJ')tz + A.{k)Ty where A(A:) = AqA: with 
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Ao a constant. The c?-vector then has just two compo- 
nents, dy = A{k),dz = ek — To evaluate W for this 
system, note that A{k) for the Hamiltonian is given by 
A{k) = e/j — /i + iA(fc). It therefore foUows that the 
winding number associated with the phase of the quan- 
tity, z{k) = Det{A{k))/\Det{A{k))\ = exp(i6'(fc)) is the 
same as the winding number N of the phase of the unit 
vector d{k) on the d — 1 Brillouin zone. The invariant 
W vanishes for /i < and is equal to 1 for /i > []. The 
point fi — marks a topological quantum transition at 
which W changes by 1. Note, however, that it is also 
possible for W to jump by 2 or any other integer (for ex- 
ample, when A(fc) ~ Aq sin_fc + Ai sin2A: with no relative 



phase between Aq and Ai [2l|) giving rise to a topologi- 



cal transition between two ground states with number of 
MFs at each end differing by more than one. 

Z invariant for the SO coupled nanowire: We now 
consider the case of a SO coupled semiconductor 
nanowire with a Zeeman coupling and a proximity- 
induced s-wave superconductivity [J-lg. In this case, 
from Eq. ([3]) and Eq. (O we have, Hq = (efc — /i) 
af{k)ay + Vzcfx and A — iA^ay so that A{k) = {ck — 
n) + af{k)ay + VzO'x + iA^ay. Here we have generalized 
the SO coupling term to have a general wave- vector de- 
pendence with the constraint f{k — ^ ±7r) — ^ 0. We find 
that 

Det{A{k)) = {eu - fif + A^ - V| - a^f{k) + 2iAoaf{k) 

(14) 

has a non-trivial winding number whenever the Pfaffian 
of the Hamiltonian (at fc = 0) — fi)"^ + Aq — F| — 
Q^^/^(^) < 0- To see why this is so consider a simple 
model for the semiconductor bandstructure ek=o = 
and ek=Tr/a ^ Mj^j^^- Similarly, because the Rashba 
SO coupling is a result of broken inversion symmetry 
f{k) = for fc = 0,7r. Therefore we assume that f{k) 
increases from to a maximum value and then decreases 
as k changes from to tt. In the case when the Pfafhan 
changes sign in going from fc = to A; = tt, we see that 
Det{A{k)) changes from being a negative real number 
(i.e. with argument Arg{Det{A{k — 0))) = tt) to a pos- 
itive real number (i.e. with argument Arg{Det{A{k — 
tt))) = 0) taking a route in the complex plane over the 
positive real axis (i.e. Im{Det{A{k))) > 0). Thus, from 
an inspection of this trajectory, the winding of the phase 
of Det{A{k)) i.e. Arg{Det{A{k))) must equal tt and the 
winding number W = 1. Since the nanowire is in the 
class BDI, the number of MFs at each end (and the in- 
variant W) can be any positive inte ger (i.e, not just or 



Z2 invariant is trivial in this case and is not adequate to 
describe the topological properties of the nanowire. 

Z invariant for the two-channel SO coupled nanowire: 
Consider a a quasi-lD semiconductor wire with two pairs 
of relevant confinement- induced bands (two channels). 
The BdG Hamiltonian for this system proximity-coupled 
to an s-wave superconductor is given by the appropriate 
generalization of Eq. [3l written as 

H2ik) = (ek - (/i + 6fi) - Sfipz)T;, + Vza-xTz + akayT^ 
+ (Ao + Ai2Px)cryTy, (15) 

where Px.y.z represents Pauli matrices that account for 
the inter-band degree of freedom. The confinement en- 
ergy splitting between the two bands is given by while 
the pairing between the bands is given by A 12. For this 
Hamiltonian, A{k) = (e^ — (/i -I- (5/i) — 5ppz) + af{k)ay + 
VzCTx + *(Ao + Ai2Px)'^y Let's first consider no inter- 
channel pairing potential, A12 ~ 0. In this case, the ma- 
trix A{k) commutes with pz, so that Arg(Det(A{k))) — 
Arg{Det{Ap^^+i{k))) + Arg{Det{Ap^=^i{k))), where 
_ Ap^=p{k) are 2x2 matrices similar to the single-channel 
case with the channel index fixed to = p = ±1. The 
winding number of Det(A{k)) is just the sum of the wind- 
ing numbers for each of the two channels pz — p — ±1, 



1). Similar to the p-wave case we can con- 

struct an example of such a Hamiltonian (with W — 2) 
by including the next-nearest-neighbor terms in the dis- 
persion and the SO coupling. Below we discuss a more 
realistic situation characterized by 2 MFs at each end of 
a quantum wire (W = 2) the stability of which is pro- 
tected by the chirality symmetry. The commonly used 



which is given by ^ [sgnV^j? 



^0 



{H + Sii + Sppz)^ + 1]. 



For \Vz\ > min[v'A2 -h^2^ ^A^ + {p - 2Sp)^] the total 
winding number can become 2 so that there are two MFs 
at end of the wire. This occurs for Zeeman splittings 
larger than 2|(5/i|, Vz > p, p — 5p. 

We see above that a pair of MFs at each end is ex- 
pected when the channels are decoupled or as has been re- 
cently shown at non-generic phase-transition points. 31 1 
However, a non-zero value of the topological invariant W 
clearly shows that the pair of MFs at each end remain at 
exactly zero energy even for non-zero inter-band pairing 
A12 as long as A12 is not large enough to close the bulk 
superconducting gap. Although intuitively clear, we have 
checked this numerically by computing the invariant us- 
ing Eq. [8] The Majorana nanowire in this case has two 
MFs at each end of the wire confirming the validity of 
the topological class BDI. 

Conclusion: We show that the Hamiltonian of a SO 
coupled topological semiconductor nanowire is in the 
topological class BDI with an integer Z topological in- 
variant. The Z invariant counts the number of possible 
zero-energy MFs at each end of the nanowire. We show 
that the familiar Z2 invariant of this system only gives 
the parity of the integer invariant. In contrast to the 
prevalent view, we find that under appropriate condi- 
tions multiple MFs (equal to the Z invariant) at each 
end of the wire can all be at zero energy protected by the 
chirality (reality) symmetry of the Hamiltonian. 
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